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Fritsch  and  Carlson  [3]  doTslopad  an  algoritba  wkieb  prodaees  a  aoaotoiie 
piseavise  eabie  iatarpolaat  to  a  ■onotoae  faaetioa.  Ws  show  that  the 
algorltla  yields  a  third-order  approxiaatioa,  while  a  aodifieatioa  is 
foarth-order  aeoarate. 
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1*  latrodMtioa. 

In  aAditioa  to  boi&g  a  good  approzlaatioa  to  tba  faaetioa.  it  is 
oftOB  dosirable  that  an  intorpolaat  roprodnoas  aneh  properties  as 
aonaegatiritp,  Boaotonieity,  aad  ooaTezity.  la  this  paper*  we  analyse 
three  algorithas  which  prodaoe  aoaotoae  pieoewise  enbie  iaterpolaats  to 
a  aoaotoae  fnaetioa. 

Siaee  the  iaterpolaat  is  a  pieoewise  oabie*  one  woald  hops  that  saeh 
an  algoritha  would  yield  a  third-’  or  fourth-order  approziaatioa  wheaewer 
the  fnaetioa  interpolated  is  suffieieatly  saooth.  Howewer*  if  the 
algoritha  (considered  as  a  aap  froa  the  set  of  aoaotoae  fnaetioas  to  the 
set  of  aoaotoae  pieoewise  eubies)  is  linear*  then  it  is  at  best  first- 
order  accurate  (see  de  Boor  aad  Swarts  [2]).  Coasequeatly*  if  greater 
accuracy  is  desired,  the  algoritha  aust  be  aoaliaear. 

Fritsch  aad  Carlsoa  [3]  proposed  such  an  algoritha.  Girea  an  initial 
pieoewise  cubic  iaterpolaat*  they  aodify  the  deriTatiwe  walues  of  that 
iaterpolaat  (where  aecessary)  to  prodaoe  a  aoaotoae  pieoewise  cubic 
iaterpolaat.  Siaee  the  aodifieatioa  process  is  aoaliaear,  one  night  hope 
that  the  Fritsch-Carlsoa  Algoritha  is  aore  than  first-order  aecurate. 

la  Sectioa  2,  we  rewiew  the  Fritsch-Carlsoa  Algoritha  and  present  two 
aodifieatieas*  the  Two-Sweep  aad  Ezteaded  Two-Sweep  Algorithas,  lAioh  also 
prodaoe  aoaotoae  pieoewise  eabic  iaterpolaats.  la  Section  3*  we  prows 
that  all  three  algorithas  yield  third-order  approziaatioas  to  a 
aoaotoae  fanetioa.  lowewer*  ia  Seetioa  4*  we  deaoastrate  that  neither  the 
Fritsch-Carlsoa  Algoritha  aer  the  Two-Sweep  Algoritha  is  a  fourth-order 
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Mtliod*  wh»xc»  is  the  eese  of  the  letter  elgoritha,  we  seenie  that  the 
iaitial  approxiaate  derlTatiwe  walaea  are  set  fonrth-order  aeearate.  On 
the  other  hand*  the  Extended  Two-Sweep  hlgoritka  ic  a  fonrth-order  aethod 
-if  the  initial  approxiaate  deriwatiwe  walaea  are  third-order  aeearate. 
Finally,  aeae  aaaerieal  exaaples  are  preaented  in  Section  5. 

For  brevity  and  eiaplieity,  we  asanae  that  the  fonetion  interpolated 
is  aonotone  inereasing  thronghont  the  reaainder  of  the  paper.  The 
extension  to  deoressing  fnnetions  is  trivial. 

2«  Algerithas. 

Xn  this  section,  we  review  the  Fritsehr-Carlsoa  Algoritha  and  present 
two  aodifieations,  the  Two-Sweep  and  Extended  Two-Sweep  Algorithas. 

The  basis  of  the  Fritsch-Carlsoa  Algoritba  is  a  techniqne  for 
deteraiaiag  whether  a  onbic  polynoaial  p(x)  is  aonotone  on  the  interval 
.  Central  to  this  techniqne  is  the  closed  region  M  (see 
Fignre  2-1^)  bonnded  by  the  axes  and  the  *qpper  half'  of  the  ellipse 

X*  ♦  y*  4-  xy  -  6x  -  6y  +  9  ■  0.  (2.1) 


Also  shown  in  Fignre  2-1  are  the  closed  regions  A, ...,£  nsed  in  the 
ei^ression  of  the  algerithas.  A  segaent  of  the  line  x  4  y  >  4  foras  the 
border  between  the  regions  A  8  *>4  between  the  regions  and  £. 
The  region  £  is  bonnded  by  the  lines  x  ■  S  and  y  ■  9. 
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Frit»eh  aad  Caxlson  [3]  tlio*  tlist  pCz)  is  boboIom  ob  0BI7 

if  (p*(Zj)»p'(Xj^j^))  s  Bbsrs 

*  K*3^  "  (  (zA^tPsi^)  :  (ztp)  s  K  ), 

Ai  -  -  p(*j)l/kj.  -  z^^j  -  Zj. 


FigBZS  I"!:  Tks  aoBOtoBicity  rsgioB  K  bbA  Bssoeistsd  eztszioz 
zsgiOBS  A* ••••£•  kll  zBfioBS  zre  elossd. 


Is  also  seals  tks  rsgieas  bp  A^  aad  rsfsr  to  tbsa  as  Aj •••••&£> 

rsspsetiTslp.  Bowstsz,  if  A^  •»  0»  ws  sztsad  tbis  eoBTSBtiOB  bp  t»iBg 
to  bs  tbs  wbols  first  qaadram;  all  otbsr  rsgioBs  eoBtraet  to  sitbsr  poiats 
or  liass  ia  tbs  obrioas  wap. 
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HkM*  startlag  witk  a  aat  of  faaetioa  Talaas  {f(z^))  aad  approziaata 
dariTatiTC  Talaaa  {4^)«  it  la  aaap  to  datazaijw  vkatkar  tha  eabic  Hazalta 
iatarpolaat  of  tbaaa  Talaas  la  BoaotOBa.  MoraoTaz*  if  tba  iaitial 
iatarpolaat  is  aot  moaotomt,  thaa  tfca  aeaditioa  oa  p*  iadieatas  bar  tba 
Talaas  {d^}  sboald  ba  aodifiad  to  aaka  it  aoaotoaa. 

3 

Figara  2-2  prasaats  a  tbraa  stap  aata-algoritba  for  fiadiag  a 
aoaotoaa  cable  Haraita  iatarpolaat.  Oalp  Stap  2  is  spaaifiad  eoaplatalp. 
la  Stap  1*  aay  taebai^aa  for  eoapetiag  tba  iaitial  approziaata  darlTatlTa 
Talaas  (d^}  is  aoeaptabla.  altboagb  tba  aoearaep  of  tba  iaitial  Talaas  is 
oaa  of  tba  priaa  factors  ia  dataraiaiag  tba  aeoaracp  of  tba  iatarpolaat. 
Tbraa  possible  iaplaaeatatioas  of  Stap  S  arc  daralopad  ia  tba  raaaiadar  of 
this  saetioa. 

Stap  1:  Coapata  tba  iaitial  approziaata  dariTatiTa  Talaas  (d^). 

Stap  2:  Easara  that  aaeb  d^  is  aoaaagati^a. 

FOE  i  :*  1  81EP  1  WIXL  a  DO 
d^  aaz(d^»0); 

Stap  3:  Hodifp  (d^)  so  that  aaA  ordered  pair  •  K^* 

Figara  2"S:  Prallaiaarp  Algoritba. 


Altboagb  Stops  2  aad  3  eaa  ba  ao^iaad  aasilp  aaTiag  oaa  pass  tbroagb 
tba  data,  aoasidariag  tbasa  two  steps  aaparatalp  siaplifias  tba  aaalpsis. 


If  Step  3  tezmiutes*  thea  tSe  algoritta  prodaeet  a  set  of 
appzoziaate  deriTatlTo  aaluee  vhiek*  togetket  vitk  tke  faaotioa  Talaas 
(f(z^)].  detexaiae  a  aoaoteae  eable  Hezaite  latezpolaat  of  f.  Tke 
difflealty  ia  iapleaeatiag  Step  3  ie  tkat  aodifyiag  oae  deziratiTe  aalae  d^ 
affeeta  both  of  tke  ozdezed  paiza  Beoaaae  of  tke 

akape  of  X*  deezeaaiag  tke  aagaitade  of  d^  la  aoriag 
fozee  ^i-1*  Tioe  Tozaa. 

Poz  tkia  reaaoa,  Fzitaek  aad  Cazlaoa  baae  tkelz  algozltba  oa  a  zegioa 
X  propezly  eoataiaed  ia  X  vfth  tke  folloviag  i^oztaat  pzopezty: 

If  (x,y)  a  X  aad  0  i  x*  i  *  ^  L  7*  i  7*  tkea  (x'*^,y'*^ )  a  X* 

Tke  Fzitaok-Cazlaoa  Algozltba  eoealata  of  Stepa  1  aad  2  of  tke  Pzelialaazy 
Algozltba  togethez  with  Step  3  aa  akowa  la  Flgaze  2-3.^ 

Altezaatlwely,  aay  teekalqae  foz  pzojeetlag  the  polata  (^£'^£4.2^ 

X£  wkleh  la  gaazaateed  to  tezalaate  eoald  be  aaed  la  Step  3.  Oae  each 
aetkod,  the  Two-Sweep  Algozltba,  la  akowa  la  Flgaze  2-4. 

Oa  the  Fozwazd  Sweep,  oaly  tke  aeooad  ooapoaeat  of  eaeb  ozdezed  palz 
la  altezed,  ao  tkat  aodlfylag  ^^£*4^4.2)  4oea  aot  affect  foz  j<l. 

Coaaeqaeatly,  it  la  eaay  to  aee  tkat  a  X£U  &£  u  §£  aftez  the 

^  Beze,  agala,  we  kawe  aaed  tke  aotatloa  X£  to  ataad  foz  X*A£. 
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St«p  3:  Modify  (d^)  so  thst  oseli  ordarod  pair 
FOR  i  :-  1  STEP  1  IIRTL  a-1  DO 
IF  («£.«£+£)  i  i£  XHBJ 

Coapots  d^  sad  d^^^^  so  tMst 
(s)  0  i  dj  i  d^. 

(b)  0  i  dj^j  i  d^^j,  sad 

(«)  (*£.*£+1)  • 

•“  'i+1  *“  Cr 

Figaro  2-3:  Step  3  of  tMe  Fritsefc^Csslsoa  Algoritta. 

Step  3:  Modify  (d^)  so  thst  esoh  ordered  pair  (d^rd^^^)  s 

Forward  Sweep  -  aodify  the  seeoad  ooapoaeat  oaly. 

FOR  i  :•  1  SIEP  1  UNTIL  a>l  DO 

IF  <«£.«£+£)  «  Cj  raw 

d  :*  3A  2 
®i+i  “i* 

ELSE  IF  (4£.S£4.£}  •  A£  u  &£  IBEN 

Decrease  to  project  t^£»A£4.£)  oato  the  boaadary  of  K£ 

Bsckward  Sweep  -  aodify  the  first  eoapoaeat  oaly. 

FOR  i  a>l  STEP  -1  XMTIL  1  DO 
IF  (<£*3^4.£)  •&£»&£  THEN 

Decresso  d^  to  project  (^£*^£4.£)  oato  the  boaadary  of  ](£* 


Figare  2-4t  Step  3  of  the  Two-Sweep  Algoritha 


Forward  Sweej. 
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On  the  Backward  Sweep,  only  the  first  eoaponent  of  each  ordered  pair 
is  altered,  so  that  aodifyiag  does  not  affect 

Horeower.  decreasing  the  nagnitnde  of  d^  ensnres  that  the  neighboring  point 

(dj_j^,d^)  reaains  in  o  |i^_j  y  ^-1*  **  •“ 

into  K^.]^  ^7  decreasing  the  nagnitnde  of  d^_^^  on  the  next  pass  throngh  the 

loop.  Therefore,  after  the  Backward  Sweep  is  eonpleted, 

the  associated  cnbic  Hemite  interpolant  is  nonotone. 

The  najor  short*-eoaing  of  the  Two-Sweep  Algorithn  is  that  it  nay  nowe 
a  point  anoh  farther  than  necessary  when  projecting  it  into 

This  problen  is  nost  acnte  in  the  regions  A  ^d  £  close  to  the  points  (0.3) 
and  (3.0).  respeetiwely.  where  the  honndary  of  M  is  tangent  to  the  axes 
(see  Section  4).  Therefore,  we  now  consider  the  Extended  Two-Sweep 
Algorithn  described  in  Figure  2-5. 

If  the  ordered  pair  of  approxinate  deriwatiwe  walnes 
not  lie  in  then  this  algorithn  allows  the  nagnitnde  of  d^  to  be 
increased  on  the  Forward  Sweep  and  the  nagnitnde  of  d^^^  to  be  increased  on 
the  Backward  Sweep.  Howerer.  the  anonnt  by  which  they  can  be  increased  is 
constrained  by  the  reqnirenent  that,  on  the  Forward  Sweep,  the  preceding 
ordered  pair  renain  in  y  y  and.  on  the 

Backward  Sweep.  ***^  renain  in  Beeanse  of  these 

constraints,  it  is  clear  that  after  the  two  sweeps  of  the 

extended  algorithn  have  been  eonpleted.  Consequently,  the  associated  cubic 
Bemite  interpolant  is  nonotons. 
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St«p  3 
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:  Modify  {d^}  to  that  oteli  ordered  pair  *  K^* 

Forward  Sweep  -  aodify  the  aeeoad  eoapoaeat  oaly  aaleaa  the 
ordered  pair  liea  ia 

FOB  i  1  STEP  1  WTIL  a-1  DO 

CASE  (d^.d^^j)  a  C^: 

^i+1  3^i* 

CASE  (d^,d^^^)  a 

Deereaae  d^^^^  to  project  oato  the  bonadary  of  Mj* 

CASE  (dj.d£+i)  t  Aj.* 

laereaae  d^  aatil  either 

(a)  reaohea  the  bonadary  of  or 

<h)  reaehea  the  bonadary  of  u  u  ®i-l 

(if  i  >  1); 

IF  i  Mj  IHEN 

Deereaae  d^^^  to  project  tdj»dj^j^)  oato  the  bonadary  of 

Backward  Sweep  -  aodify  the  firat  eoapoaeat  oaly  aaleaa  the 
ordered  pair  liea  ia 

FOB  i  :•  a-1  STEP  -1  1MTIL  1  DO 

CASE  (djrdj^j)  a  I)^: 

Deereaae  d^  to  project  (d^.d^^j^)  oato  the  bonadary  of 
CASE  (dj.dj^j)  a  E^: 

laereaae  d^_^^  aatil  either 

(a)  ^d^,d^^^)  reaohea  the  bonadary  of  or 

(b)  Ui+i»*i+2^  reaehea  the  bonadary  of  Ej+j  (if  i  <  a-1); 
IF  (d^.dj^j)  4  THEN 

Deereaae  d^  to  project  (d^^d^^^)  oato  the  bonadary  of 


Figara  2—5:  Step  3  of  the  Ezteaded  Two-Sweep  Algorithm 


3.  mrtf-Ord«r  CoBTcrgcufc 
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In  tlilt  section,  ve  pzoTS  tlist  sseh  of  the  slforlttas  presented  in 
Section  2  yields  s  third-order  spproxiastion  to  s  aonotone  fnnetion, 
provided  that  the  initial  approziaste  deriTative  valnes  are  seoond— order 
aoenrate  and.  in  the  ease  of  the  Fritseh-Carlson  Alporitha.  that  £  is 
snitably  chosen. 

Ye  hepin  by  considering  vhat  restrictions  on  the  region  £  are 
necessary  for  the  Fritseh-Carlson  Algoritha  to  be  third-order  aoenrate.  To 
this  end.  the  following  result  is  nsefnl. 

Lsaaa  3^1:  If  p^(z)  and  P2(x)  are  two  pol3raoaials  of  degree  three  or 
less  that  satisfy 

PjCzj)  -  pjCXi)  and  » 

then 

aaz  {  IpjCx)  -  PjCx)!  :  i  x  i  Xj^2  )  (3.1) 

^  6/1  ^  "  Pi(*i)l.  Ipi<*i+1^  "  '2^*i+l^' 


Proof:  Ewalnating 


.x-x.  .,,2 


Pl^x)  -  P2<*>  “  (x-Xj)  ["  “  ®2^*i^^ 

*  ^*“*i+i^  [  h^^l  ^*i**i+i^  "  ®2^*i+i^^ 

at  the  points 


^i  -  *1  ♦  ^2  ’  *1  -  *1  ♦  I5  ♦ 


-  10  - 


yield* 

»l<Tj)  -  Pj(Ti>  -  4^  t  I5  •*  Tn*  ■  »2‘*l” 

-  ‘I  -  ylsr'  > 

•nd 

^2  ^  JH ^ 

respeetlTely.  If 

Ip^Cxj)  -  P2<*i)l  i  •Pi<*i+i>  “  P2^»i+l^l* 

than*  froa  (3.2), 

IpjCyf)  -  ^  ^  iP2<*l>  "  P2^*i^* 

"  ^2  "  7n^  "  **2**i+l^*  ^ 

2.^  lpi<V  -»£<»j>l. 

vhieh  iaplies  (3.1).  On  the  other  head,  if 

Ipi(Xi)  -  P2^*i^l  i  ^^'i^^i+l^  "  **2^*i+l^*' 
then  (3.1)  follows  froa  (3.3). 


(3.2) 


(3.3) 


Dales*  (1,1)  s  2  (the  elosare  of  fi)*  t^e  FritsehrCerlsoa  Algoritha  is 
St  best  first-order  soearste.  Coasidor  the  spproxiastloa  to  f(z)  ■  z  oa  s 
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•aifexa  ■•sk.  la  thit  ease. 

f'(Xj)  «  -  Aj  -  1,  for  1*1,..., 8-1. 

Coasoqaoatly.  fox  axoA  i.  oaa  of  tad  d^^^  net  be  boaadod  away  froa  1. 
aad  tba  raaalt  follows  froa  Loaaa  3.1. 

8i>llaaly*  aalass  X  c  £*  wbara  X  Is  tbe  eloaed  triaagle  with  wertieas 
(0.0).  (2.0).  (0.2).  the  Pritaeh-Caxlaoa  Algoritha  it  at  beat  aaeoad— order 
aeeazate.  Assaae  soae  poiat  (a.2-s).  0  i  a  <  1.  oa  the  'i^per  half'  of  the 
hypoteaese  of  X  1*  aot  ia  £  aad  ooaaider  the  approxlaatioa  to  f(x)  ■  (x-a)^ 
oa  the  iaterwal  la.bl.  For  aay  h  i  a  (b-a) .  ehooae  a  set  of  kaots 

(Xj)  aad  an  iateger  J  saeh  that  Xj  •  a  +  hj  •  h  •  aax{h^).  With 

this  ehoiee  of  Xj  aad  h^.  x^^j  i  b. 

£'(x.)  2(x.-a)  ^'<*4+1 >  3(x,-a)+2h, 

-ip  ■  Jupirs;  ■  -  - 

lloreoTor*  A^  2  Therefore,  whea  the  Fritsch-C^rlsoa  Algoritha 

teraiaates.  at  least  oae  of  the  approxiaate  deriwatiwe  walaes  (d^)  aast 
satisfy 

If'(Xj)  -  d^l  I  eh. 

for  soae  eoastaat  e  >  0.  aad.  by  Lsoaa  3.1.  the  aasoeiated  eabie  Heraite 
iaterpolaat  is  at  best  seeoad-order  aeearate.  A  aiailar  resalt  holds  for 
the  'lower  half'  of  the  hypoteaese  of  X> 


Ob  the  other  head,  if  X  c  ft»  thea  the  Fritseh-Caxlaoa  Algoritha  is 
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thixd-ordar  aeoaxate.^  Bafoxa  pxoriag  this  xaaalt,  wa  atata  tha  folloviag 
aaaful  laaaa. 

LumB  3'»2:  If  f  a  C^Ia*b]  la  Boaotoae  iaexaasiag.  than*  fox  aap  of 
tfca  alfoxittea  daaexibad  la  Saetloa  2, 

d*  i  0  aad  If'(Xj)  -  d*\  i  lf»(Xj)  -  d^l,  1-1 . a, 

wbaxa  d^  aad  d^«  xaspaetlTaly*  axa  tba  appxoziaata  daxlvatlTa  Talaat  bafoxa 
aad  aftax  tba  azaeatloa  of  Stop  2. 

Proof:  If  d^  la  aodlflad  la  Stop  2,  tbaa  d^  <  0  aad  d^  -  0  (aaa 
Flgaxa  2-2).  Haaoa,  alaea  f'(z^)  2.  0, 

If'(Zj)  -  dji  -  lf'(Zj)l  <  lf»(Zj)  -  djL 

Oa  tba  otbax  haad»  If  d^  la  aot  aodlflad*  tbaa  1  0.  ttvS.!). 

Ikooraa  S*S:  Aaaaaa  that 

1.  f  a  C^[a,bl  la  aoaotoaa  laoxaaalagi 

2.  tba  laltial  daxlvatlra  appxozlaatloaa  {d^}  aatlafp 

If'(Zj)  -  dji  i  Ob*.  1-1*. ...a* 

^  Tba  foax  xagloaa  jL*...*^  eoaaldaxad  la  [3]  all  eoatala  tba 
txiaagla  I. 
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fox  SOB*  eoBttaat  e; 

3  •  Z  c  £;  aad 

4.  iAomtox  a  point  pxojoetod  into  £^>  the  new  point 

aatiafies 

2Ai  i  dj  +  dj^^, 

i.e.*  the  point  ie  not  moved  *aBeh  farther'  than  neeeeeary. 

Than  the  nodified  approziaate  derivative  valnee  (d*)  prodnoed  bp  the 
Frlteeh-Carleon  Alforithn  eatiefy 

If'(Zj)  -  d*l  i  [e  ♦  ||f^*’ljh^.  i*l,...,a.  (3.4) 

Coneeqaentlp*  the  aaeoeiated  aonotone  enbie  Beraite  intexpolant  ie  a  third- 
order  approziaation  to  f. 

Broof:  Froa  Lomb  3.2. 

d^  i  0  and  If'(z^)  -  dji  i  oh*  (3.5) 

at  the  taraination  of  Step  2. 

Aecnae  that  d^  ie  modified  in  Step  3  when  projected  to 

(dj_j,dj)  e  *i-i»  ^1-1'  *i  differ  froa  the 

^  dT  and  d?  ^  are  approziaate  derivative  valnee  that  have  been  aodified 
either  race  or^tviee.  reepeetivelp. 


initial  Talaaa  and  d^  aatiafying  (3.5)*  bat 

0  i  i  di.i  i  nnd  0  i  dj  i  d^. 

If  f*(x^)  i  d*.  than 


0  i  d*  -  £*{x^)  i  -  f'(Xj)  i  ok*. 


Tkarafoxa.  aaaana  tkat  f'(x^)  1  d^.  Note  tkat 
“i-i  - 


for  aona 


Tj.j  e  J'***  Aaan^tion  4, 


*^1-1  i  ^i-1  *  *i' 


ao  tkat 


f'(x^)  «  4 j  i 
Tkaxafora* 

0  i  f»<Xj)  -  d* 

i  Vl  ” 

i  «i_l  -  ^’<*1-1^  ♦ 

He  *  ||f^*^ljk* 

by  Aaanaption  2. 

If  d^  ia  daexaaaad  to  d^  to  prejaet  ^4^*4^^^^)  into  on  tha  naxt 
paaa  thxengk  tka  loop,  tkan  a  aiailax  axgnnaat  akowa  tkat  inaqnality  (3.4) 


xanaina  Talid 
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Esteatially  tfc*  saa«  argncat  skows  tkat  tka  Tvo-Svaap  Algoritka  it 
tklrd-order  aeoaxata.  Bevavar*  tka  Bataadad  Tv»-8*aap  Algor itka  aay 
iaeraaaa  aoaa  approziaata  dariaatiTa  valBoa.  Tkarafora*  «a  adopt  a 
diffaraat  approaek  baaad  ^oa  tkc  follovlag  lasa.^ 


Laaaa  S'«4:  Aaaaaa  tkat 

1.  f  a  C^[a,b]  la  aoaotw^  iaoraaaiag;  aad, 

2.  for  aoaa  o  >  0.  (f'<z^,j) ,f»(z^))  d  lJ_j  -  wkara  J*  la  tka 

eloaad  trlaagla  vltk  Tartleaa  (0,0) »  (2^»0)  •  (0,2-*«). 

Tkaa 


‘1-1 


(3) 


Ik 


1-1 


aad 


♦  !•(«,) 


1-1’ 


Proof:  If  <f'(xj_j),f»(x^))  4  Uaa 

(2+o)Aj_^  <  I'txj.j)  ♦  f'(xj). 


7 

la  paaalag*  aota  tkat  tkla  laoaa  eaa  alao  ba  aaad  to  prora  a  diffaraat 
▼araloa  of  Tkaoraa  3.3:  If  Aaaaaptloaa  3  aad  4  ara  raplaead  by 

3.  £  for  aoaa  a  >  0, 

tkaa  tka  Frltaek-Carlaoa  Algorltka  la  atlll  tklrd-ordar  aeearata. 
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I 

I 

\ 


■  I 


lav«T«r» 

for  •ea*  t  •*  fkat 

*1-1  <  si'‘”<^-i>‘i-i  * 

Fiaally*  Bsiag  (3.6), 

'•'•i-l>  *  ''Ol'  <  [l  *  lil*''’’'-**!-!- 

as.D. 

Ikoeroa  3*5:  Astoac  tliat 

1.  £  •  C^U.bl  it  aonotoae  iaorottiag;  tad 

2.  tbt  laltitl  dtriTttivt  tpprosiattieat  (d^)  tttitfy 

l£'(Xj)  -  d^l  i  ck^,  i»l.....a 
fox  toat  eoasttat  e. 

Tkta  tkt  aodiflod  tpproxlatte  dtriTttlTO  Ttlaot  {d*}  prodaetd  by  titktr  the 
Tvo>8voep  or  the  Exteaded  Two-Soetp  Alforitka  tttltfy 

|f»<Xj)  -  d*l  i  atxCc,  jlf^’^ljh*.  W,...»a.  (3.7) 

Coateqaoatly,  the  attoeltted  aoaotoae  oable  Horaite  iaterpoltat  it  t  third- 
order  tpproxiattioa  to  f. 

Proof t  By  Ltaat  3.2«  the  tpproxiatte  deriTttiTt  rtlaot  tttitfy 

2  0  It'(Sj)  -  i  a* 


k 
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•t  tli«  eoaplfttloa  of  Stop  2  of  oltkox  olporitte.  TSoroforo*  tLop  also 
satlafp  (3.7).  Bslov*  vs  akov  tkat*  if  all  tka  approzlaato  doriTatlTO 
▼alaoa  aatiafjr  (3.7)  vksa  om  la  aodiflod  la  Stop  3*  tkoa  tko  aodlflod 
▼alas  also  aatlafloa  (3.7).  Tkat.  tko  tkooroa  folleoa  bp  ladaetloa. 

la  tko  Eztoadod  Two^Svoop  Algorltba*  d^  la  aodlflod  la  Stop  3  oalp  If 

1.  projoetod  dovaoarda  la  tko  Porvard  Sooop. 

2.  (^£*^£4.2^  projoetod  to  tko  rlgkt  la  tko  Porvard  Svoop, 

3.  projoetod  to  tko  loft  la  tko  Baekvard  Swoop*  or 

4.  projoetod  ^vards  la  tko  Baekvard  Swoop. 

Por  tko  Tvo-Svosp  Al|oritha«  oalp  Caaoa  1  aad  3  are  applieablo.  Tkoroforo* 
proTlag  (3.7)  for  tko  Eztoadod  Tvo-Svoop  Algoritka  also  akeva  tkat  tkla 
laoqaalltp  la  valid  for  tko  Two-Swoop  Algoritka. 

Coaaldor  Caao  1  first:  Pi^ojoetod  dovawards  la  tko 

Porvard  Swoop.  If  f'(z^)  ^  d^*  tkoa 

0  i  d*  -  f'(Zj)  i  Sj  -  |lf^*^ljk*, 

alaoo  d^  X  Tkoroforo,  aaaaao  tkat  f*(z^)  2. 

(f'(*j_j)  .f’(x^))  a  tkoa 

f'(Zj)  i  SAj.j  i  dj. 

a  eoatradietioa.  Tkas,  (f '(z^_^) ,f'(z^))  i  lAoaeo 

f'(Zi_i)  ♦  f»(Zj)  < 


by  3.4.  Siaee  f(*^)  i.  d*  1  0  aai  both  aad  £'(x^)  are 

aoaaogatiT*. 

0  i  f'(Xj)  -  d*  i  £»(xj)  i  f*<xj_j)  ♦  f'(x^)  i  |lf**^l,h*. 

Noxt  eoxsidor  Case  2:  yrojeeted  to  the  right  is  the 

Forward  Sweep.  If  d^  i  f*(x^).  then 

0  i  f(Xj)  -  d*  i  f'(x^)  -  dj  i  »ax(e.  |lf**^l,)h*. 

aixee  d^  ^  ^i*  Therefore,  asawae  that  d^  2.  f*(x^).  If 
(f»(x^).f*(xj^j))  4  ij'*.  thea 

by  IrOBBa  3.4,  Bat  d*  i  jA^  alaae  a  »®  that 

0  14*-  f(,,)  1  4*  1  |4,  1 

Ob  the  other  head,  if  (f'(x^) ,f»(Xj^j))  a  thea 

f'(Xj)  +  f'(xj^j)  i  jAj  i  dj^j  -  jAj, 

aiaee  a  4^  iapliea  that  3A^  ^  ^i4>l*  Ae-srraagiag  teraa. 

|4,  ♦  1  4.^1  -  f 

i^aee 

df  +  f'(Xj)  i  dj^j^  -  f»(x^^j), 

aiaee  d^  ^  |a^.  Therefore, 


Caaet  3  and  4  are  haadlad  ia  a  alailax  aaaaax 


O.'B.O 


4.  Faartk-Qrdax  Comnx$mmom^ 

la  thia  aaetioa.  «a  daaoaatxata  tkat  aaitkax  tka  Fxltaek-Caxlaoa 
Algoxltka  aox  tka  T«o-S*aap  Algoxitka  la  a  foaxtk-oxdax  aatkod.  vkaxa,  ia 
tka  oaaa  of  tka  lattax  algoxitka.  va  aaaaaa  tkat  tka  iaitial  appxoxiaata 
daxiTatiaa  aalaaa  axa  laaa  tkaa  foextk-oxdax  aecaxata.  Oa  tka  otkax  kaad. 
tka  Extaadad  Two-Svaap  Algoxitka  ia  a  foaxtk-oxdax  aatkod  if  tka  iaitial 
appxoxiaata  daxiratiTa  aalaaa  axa  tkixd-oxdax  aecaxata. 

To  aaa  tkat  tka  Fxitaek-Caxlaoa  Algoxitka  ia  aot  a  foaxtk-oxdax 
aatkod.  eoaaidax  tka  faactioa  f(x)  •  (x-1)  oa  tka  lataxral  [0.3] .  Fox  aay 
poaitiTa  iatagax  a.  let  tka  kaota  ka 

Xj  ■  3ik.  i  -  0*1.....3a+2.  akaxa  k  -  . 

A  aiapla  eoapatatioa  akowa  tkat 

-  f**-  ''‘vi’  ■  “*  *.  •  h*' 

akaaea 


(4.1) 


ia  oa  tka  koaadaxj  of  Oa  tka  otkax  kaad.  aap  xagioa  aaad  ia  Stap  3  of 
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the  Friteeh-CarleoB  Algorltha  aost  be  eoatelned  1b  the  regioB  the 

eqaare  vith  Tertioes  (0,0)  (0,3),  (3,3),  (3,0),  eo  that  the  aodified 

•  0 

derlTBtiTe  approziaetioB  d  aost  aatiafy  d  ^  3A  .  Thaa, 

■  B  B 

i  f(.,)  -  >A,  -  ik*. 

aad,  fxcMB  Laaaa  3.1,  the  Fritaeh-CazlaoB  Algoritha  yields  at  best  a  third- 
order  approziaatioB  to  f. 

To  see  that  the  T^o-Sveep  Algoritha  is  aot  a  foortb-order  aethod  if 
the  laitial  approziaate  deriTstiTe  Talaes  are  less  thaa  foarth-order 
aoearate,  oaee  agala  eoasider  the  faaetioa  f(z)  >  (z-1)  oa  the  iaterral 
[0,31.  For  2  1  p  1  4,  ehooae  the  kaots  (z^)  aaeh  that,  for  soae  j, 

f»(Zj)  -  3h**,  •  3lh*  -  22^"*^'*  +  h®), 

aad 

Aj  -  h*  -  3h^*®^*  +  3h*. 

It  is  easy  to  cheok  that  (f '(z^) ,f'(Sj^j^))  is  oa  the  boaadary  betweea  Hj 
aad  Aj  ud  that  (f'(z^) ,f'(z^^^))  a  3or  i  #  J.  Let  dj  •  0  aad 

d^  ■  f'(z^)  for  i  d  J.  Thaa  dj  is  a  p^^-ordet  approziaatioa  to  f'(Zj)  aad 
all  other  d^  are  ezaet.  la  additioa,  aiaee  dj  <  f'(Zj),  it  follows  that 
(dj.dj^^j)  a  Aj'ij  “4  **l'^i+l^  •  ^*^1  J*  Coaaeqaeatly,  the  oaly 

approziaate  deriratiwe  wslae  that  is  aodified  by  the  Two-Sweep  Algoritha  is 
dj^2  aad  it  is  set  to  ■  3A^  oa  the  Forward  Sweep.  Heaee, 
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mnd,  by  Laaia  3.1«  tb«  Tvo-Svaap  Algoritta  yialda  at  baat  as  ordar  2*^ 
approziaatloB  to  f.  la  paztiealaz*  if  tka  Tvo^Saaap  Algoritte  ia  aaad  to 
■odify  tka  dariTatiTa  Talaaa  of  a  eobic  apllaa  latarpolaat,  tkaa  tka 
raaaltiag  aoaotoaa  piaeawiaa  eabic  iatarpolaat  aay  ba  of  ordar 
ratkar  tkaa  4,  aiaea  tka  iaitial  approziaata  darivativa  Talaaa  ara  oaly 
tklrd-ordar  aeearata. 

HowoTar.  for  both  tka  Fritaek-Carlaoa  aad  Two-Saaap  Algorithas,  tkia 

dagradatioa  ia  tke  ordar  of  tka  approziaatioa  ariaaa  oaly  aadar  Tory 

apaelal  eiraaaataaaaa.  If  tka  ragiea  £  aaaooiatad  vitk  tka  Fritaek-'Carlaoa 

Algoritha  eoataiaa  a  triaagla  for  aoaa  a  >  0,  tkaa»  aaiag  aa  argaaaat 

aiailar  to  tka  oaa  aaployad  ia  tka  proof  of  Tkaoraa  4.1,  oaa  eaa  akow  tkat 

tke  dagradatioa  ia  tka  ordar  of  aitkar  of  tkaaa  two  algoritkaa  oeears  oaly 

ia  iatarrala  iaaadiataly  adjaoaat  to  aa  iatarral  aoataiaiag  a  root  of  f'  of 

azaet  araltiplicity  two.  Moraorar,  for  tka  Two-Swaap  Algoritka,  tka 

dagradatioa  oaeara  oaly  if,  aa  k  ->  0,  tkara  ara  iafiaitaly  aaay  grida  aaek 

ooatalaiag  aa  iatarral  [z^,z^^j^]  aad  a  poiat  t  ia  tkat  iatarral  at  wkiak  f 

kaa  a  root  of  azaet  aaltiplieity  two  aad  tka  diataaea  batwaaa  t  aad  oaa  of 

2 

tka  aadpoiata  of  tka  iatarral  ia  laaa  tkaa  e^^k^  bat  graatar  tkaa  e2k^  for 
all  poaitira  eoaataata  aad  02. 

Aaotkar  poiat  aboat  all  tkraa  algoritkaa  akoald  ba  aapkaaiaad: 
whaaarar  k  ia  aaffieiaatly  aaall,  aoat  of  tka  iaitial  dariratira 
approziaatioaa  ara  aot  ekaagad  by  aay  of  tka  algoritbaa.  Tkaa,  if  tka 
iaitial  dariratira  approziaatioaa  ara  tkird-ordar  aeearata,  tkaa  tka 
iatarpolaat  prodaoad  by  aay  of  tka  algoritkaa  ia  loeally  a  foartk-ordar 
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approziaatioa  on  aost  intarvals.  Hocaovar,  if  tha  initial  intarpolant  ia  a 
enbie  aplina,  than  tbis  additional  aaootbnaaa  ia  loat  only  at  tba  lenota 
abaza  tba  daziTatiTa  Talnaa  aza  aodifiad. 

Wa  and  tbia  aaetion  vitb  a  oonTazpanea  zaanlt  foz  tba  Eztandad  Two- 
Swaap  Algozitbn. 

Ibaoraa  4;U.:  Aaanaa  that 

1.  f  a  C^Ia.bl  ia  nonotona  inezaaaing; 

2.  wbanawaz  f*(x)  ■  £**(x)  •  £^^\x)  «  f^^^{x)  ■  0,  tbaze  ia  a  6  >  0 

amob  that*  if  y  a  [x.x'fb)  n  Ca»b]>  than  aitbaz 

a.  f'(y)  -  0  oz 

b.  tbaza  axlat  oonatanta  a^,  and  z  aneb  that 

■i(r-x)*  i  fir)  i-jCr-*)*. 

•kaza  ifi^ijo*®**^** 

and,  if  y  a  (x-6,x]  n  [a,b],  than  aitbaz 

a.  f'(y)  ■  0  oz 

b.  tbaza  axiat  conatanta  n^,  n^,  and  a  aneb  that 

■3<*-y)*  i  fir)  i*4(x-y)* 

wbaza  JJ  i  ^  i  Jo  •«4  •  2  S» 

3.  tba  initial  daziratiwa  appzoxinationa  (d^)  aatiafy 


Tlieii.  for  h  tttfflclently  nail,  the  Modified  approziaate  deriTatiTe  yalBes 
{d*}  prodaoed  by  the  Extended  Two-Sveep  Algoritha  aatiafy 

If'(x^)  -  d*l  i  oh^.  l-l.....n.  (4.1) 

where 

c  -  MaxCglf^^^l^,  (4.2) 

Coneeqnently.  the  aaaociated  Monotone  onbic  Hernite  interpolant  is  a 

g 

fonrth-order  approziaation  to  f. 

Proof:  To  prove  this  resnlt.  we  eoabine  a  eoapaetness  argaaent  with 
indnetion.  The  essence  of  the  proof  is  ontlined  below;  the  details,  which 
arc  straightforward  but  tedions.  are  in  the  Appendix. 

For  each  t  c  [a.bl.  we  choose  a  &^  >  0  that  deteraines  an  open 
interval  *  (t-h^.t-^^).  where  6^  depends  npon  f  in  a  neighborhood  of 
t.  Since  (I^)  foras  an  open  covering  of  the  eoapact  interval  [a.bl.  there 
exists  a  finite  snbcovering  of  (a.bl.  Moreover,  for  h  »  aax{h^} 


g 

The  proof  of  this  resnlt  reqniros  Assnaption  2.  although  we  suspect 
that  the  theorn  reaains  valid  for  any  aonotone  c  [a.bl  function.  It  is 
also  worth  noting  that  Assnaption  2  holds  for  any  piecewise  analytic 
faaetion. 


I 
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sufficiently  sasll,  eseb  iaterrsl  c  one  of  the  interrsls  of 

the  suboovering .  The  proof  relies  hesTily  npoa  exploiting  the  loesl 
properties  of  f  on  eseh  interrsl  of  the  finite  snbooTering. 

The  sotnsl  indnetion  hypothesis  nsed  is  slightly  stronger  then  (4.1): 

1.  If  c  f'<t)  -  f»»(t)  -  0,  f^^^t)  4‘  0,  snd 

t  e 

If(Xj)  -  dji  i  +  d.ScJh^. 

2.  If  c  I^.  f'(t)  -  f(t)  -  0,  f^^Nt)  ^  0.  snd 

t  s  then 

3.  Otherwise* 

If'(Xj)  -  dji  insxU,  flf^^^ljh®. 

By  Lessis  3.2* 

dj  2  0  snd  If'(Xj)  -  dji  i  oh* 

St  the  teninstion  of  Step  2.  Consequently*  the  induotion  hypothesis  is 
sstisfied  St  the  beginning  of  Step  3.  In  the  Appendix*  we  show  thst*  if 
sll  the  spproxiaste  deriwstiwe  wslnes  sstisfy  the  hypothesis  when  one  is 


■odified  in  Step  3*  then  the  aodified  wslue  slso  sstisfies  the  hypothesis. 
Thus*  the  theoren  fellows  by  induotion.  (LCaB 


5.  Ilv*rl««l 


la  this  section,  ve  ooBpsre  tke  piecewise  eabic  iaterpolsats  prodaced 
by  COBSPL  [1] .  the  Fritscb-Csrlson  Algoiritba,  sad  tbe  Extended  Two-Sweep 
Alforitba  for  tbe  two  sets  of  aoaotoae  date  giwea  ia  Section  S  of  [3] . 

la  tbe  esse  of  CDBSPL,  we  ased  tbe  *aot-s-kaot'  boaadsry  conditions 
to  ooaplete  tbe  specification  of  tbe  eabic  spliae  iaterpolsnt.  Since 
CDBSPL  is  based  apon  s  foartb-oxder  liaesr  slgoritba,  it  does  not,  ia 
general,  prodnee  s  aoaotoae  spproxiastloB  to  s  set  of  aoaotoae  data. 

We  iapleaented  tbe  Fritscb-Csrlson  Algoritba  described  ia  [3]  sad, 
followiag  tbeir  snggestion,  we  took  tbe  region  £  reqaired  ia  Stop  3  to  be 
£2*  tit*  iaterseetion  of  tbe  disk  of  rsdins  three  centered  at  tbe  origin 
with  tbe  first  qasdrsnt.  The  resalts  ia  Sectioas  3  aad  4  sbowe  show  tbst 
tbis  aetbod  is  tbird-order,  bat  not  foartb-order,  aecarste. 

le  ased  tbe  deriTstiwe  of  tbe  cable  spliae  iaterpolsat  prodaced  by 
CCBSPL  for  tbe  initial  deriwatiwe  spprozia^tions  reqaired  ia  Step  1  of  tbe 
Extended  Two-Sweep  Algoritba.  Since  these  spproxiaste  deriwatiwe  wslnes 
are  third-order  aecarste,  tbe  aoaotoae  iaterpolsnt  prodaced  by  tbe  Extended 
Two-Sweep  Algoritba  is  foartb-order  aecarste. 

Figare  5—1  shows  tbe  iaterpolsats  prodaced  by  CDBSPL  sad  tbe  Extended 
Two-Sweep  Algoritba  for  tbe  first  data  set  (AKHIA  3)  ia  [3] .  Figare  S-2 
shows  the  iaterpolsats  prodaced  by  tbe  Fritseb-Csrlsoa  Algoritba  sad  tbe 
Extended  Two-Sweep  Algoritba  for  tbe  ssae  dsta  set.  Figares  5-3  aad  5-4 
show  tbe  iaterpolsats  generated  by  tbe  ssae  two  pairs  of  aetbods,  bat  for 


26 


the  eeeond  dmta  set  (KFN  14)  ia  [3] . 

The  interpolsBt  pzoduoed  by  CDBSPL  ie  clearly  aot  aonotoae  for  either 
data  set  and  does  not  yield  a  'Tisaally  pleasing'  approziaation  in  either 
ease. 


For  the  first  data  set,  the  interpolants  produced  by  the  Fritsch- 
Carlson  and  Extended  Two-Sveep  Algorithas  differ  significantly  on  the 
intexTal  [11,15].  Becaase  the  Extended  Two-Sveep  Algoritha  projects 
approximate  deriwatiwe  valnes  onto  the  bonndary  of  K,  it  prodnees  an 
iaterpolant  with  a  zero  slope  in  this  interral.  This  is  aot  the  ease  for 
the  Fritseh-Carlsoa  Algorithm,  sinoe  it  projects  approximate  deriwative 
▼alnes  into  the  interior  of  I|.  We  leave  the  anbjeetive  qneation  of  which 
approximation  is  'visaally  more  pleasing'  to  the  reader. 

For  the  second  data  set,  the  interpolants  produced  by  the  Fritsch^ 
Carlson  and  Extended  Two-Seeep  Algoritbms  are  virtnally  indistinguishable 
at  the  resolution  of  these  plots:  moaotonicity  imposes  a  severe  eoastrai: 
ia  this  example. 
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Figar*  S*>1:  A  plot  of  tA«  iatarpotuts  prodaoad  top  CDBSPL  (dotted  eorro) 
aad  ttoo  Bxtoadod  T*o-8voop  Algeritta  (telid  earro)  for  tho  data  aat 

AKHIA  S. 
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FifBT*  5>3t  A  plot  ot  tho  iat«rpol«at»  prodaead  bp  tb«  Fritselr'Catlsoa 
Algorltta  (dotted  earro)  aad  tbt  Bxtoadod  Twe^Booop  Algoritba  (solid 
oarro)  for  tbs  data  sot  AIHIA  S. 
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FifW*  5>Sx  A  plot  of  iatorpolosto  prodoeod  bp  CDBSFL  (dottod  emrro) 
•ad  til*  Bxtoadod  Two-8««op  Algoritta  (solid  earro)  for  tbs  data  ast 

IPN  14. 


-  SI  - 
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Iw  IS«Of  of  XkOOTM  4^^^' 

la  tKls  appoadlz.  vo  eoi^lote  tk«  proof  of  Thooroa  4.1.  To  bogia*  we 
state  sad  prows  two  asefal  IsaBss. 

Um  Sai  If  f  s  C^ts.bl  sad  f*(t)  -  f"(t)  «  0  bat  f*^Nt)  ^  0  for 
soae  t  t  tbea 

■  rtT»  ^  “•*’ 

l-3r+3T 

sad 

lf‘(z. .,)  -  A J  i  ^f^^^l.(l-T)V»  <5.2) 

‘  ^  1-3t+3t*  *  ^  ‘ 

wbere  y  •  (t-z^)/b^.  Horsower.  tbs  loeas  of  poiats 

{  (  ,  .  Sll-7}^.  ,  s  0  i  T  i  1  )  (5.3) 

1-3t+3t  1-3t+3t 

is  tbs  elliptiesl  boaadsry  of  {{. 


Proof:  laeqaslities  (5.1)  sad  (5.2)  follow  froa  tbe  Taylor  series 
ezpsasioas 

f'(zj)  -  |f^*\t)T*bJ  -  |f^^^yi)A5. 
f’(Zi+i)  “  |«***<t)(3-r)*bj  +  |f^*^yJ)(l-T)*bJ.  sad 
Aj  -  |f**\t)I(l-T)^VlAi  ♦  5jtf*^^y5>(l-T)^-T^^^y4)Y^lAj. 
for  soae  y^,  y^,  jg,  y^  •  [*j**2^2^*  wslidity  of  (5.3)  is  estsblisbed 
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•asily  froB  (2.1). 


a.s.D. 


5ia2:  Asaua  tAat 


^  *  C  [a.b]  it  BOMtOfte  iseraaaiBg; 

2.  f(t)  -  £'*(t)  -  0  bat  ^  0  fot  ao>a  t  e 

3.  *  4£« 


4.  tbe  initial  dariTatine  approxiaatioaa  aatiafy 


£'<*i)ld^  and  l£'(Xi+i)  -  i  cb' 


fox  aona  eonstant  e. 


Iban,  for  tba  oaidna  d^^^  snob  that  (d^^d^^^)  a  M^n 
If'(Xi^j)  -  d^^^l  iBaxU.  J§|f(4)|^j^3^ 


A  aiailar  xaaolt  bolda  fox  (d^.d^^j^)  a 


(5.4) 


Proof:  Tbxongboot  tbia  proof*  va  nae  iaaqnalitiaa  (5.1)  and  (5.2)  of 
5.1  vitbont  a^lielt  xafaranea. 


Coaaidax  two  eaaaa  dapaadiag  190B  wbatbar  d^^^  >  f*(x^^2). 


Caaa  1:  If  dj^j  >  f'(Xj^j).  tbaa 


0  <  dj^,  -  f'(x^^j)  i  dj^j  ~  f*(X4^,)  i  Ob  , 


ainca  d^^^  i  d^^j. 


Cate  2:  If  ^  thaa  ooatider  tvo  tobeatet  depandlBg 


a  whether 

— A  <  d 

2  A.  <  d.. 

1-3t+3t 

Cate  2.1:  If  (5.5)  it  walld.  than 

<dr 


(5.5) 


1-3t+3t 


2  “1  '  “1+1* 


(5.tf) 


tinea  the  tag^ant  of  the  eurre  (5.3)  that  font  the  bonadary  between  uid 
it  an  ineraating  fanotion  of  y  in  both  the  z  and  j  eo-ordinatet. 
Conteqaantly. 

0  i  -  *1^  A,  i 

1-3t+3t 

Cate  2.2:  If  (5.5)  it  not  walid.  than  eontidar  two  tnbeatet 
depending  ^pon  whether  T  ^  3« 


Cate  2.2.1:  If  y  >  then 


l-3T+3r 


>  1. 


Thie  boand  together  with  the  obterwatioa  that  d^  (tinea 

(^i'd^+j)  •  ahowt  that 

|a.  L  — r  A,  -  d.  i  A.  -  f'(z,)  i 


1-3t+3t 
In  addition. 


1-3t+3t 


*^i  i  Cl  i  ''<*i+l^ 


wheaoe. 
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WiS' 


0  i  -  <«  i 

Case  2.2^:  AlteraatlTely.  if  0  i  t  i  |»  tk«»  th®»«  exists  a  wiqae 


C  s  [O.yI 

.2 


m: 


A.. 


*1  “  2  i 

^  1-35+3?'* 


(5.7) 


sixes,  by  assoaptioa. 


0  i  d  i  -2* - J  A 


1-3t+3t 

aad  tbe  right  side  of  this  iaeqoality  is  a  strictly  ixereasing  fonctioa  of 
T  for  0  i  T  i  MoreoTor,  sixoe  (d^.d^^j^)  s  n 


.+  ,  . 

1-35+35*  * 


hy  (5.3).  Therefore. 


0  i  *'^*i+l^  “  *1+1 

d 

-3t+3t 

i  ♦  9(T-t)Ai. 


(5.8) 


f'(x  1)  -  Aj  ♦  2  Aj  -  2  *i 

1-3t+3t*  *  1-3T+3/  *  1-35+35* 


siaee.  for  0  £.  5  1  T  1  f* 

0  i  i  »(T-t)  . 

1-3t+3t  1-35+35 

To  boaad  3(t-5)Aj.  aote  that.  forOi5iTi5“d  f*(x^)  -  d^  i  0. 
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0  i  3t(t  -  5)Ai  i  3(t*  - 

i 


2  i  2  “4 

l-3T+3y*  *  l-S5+3r 


1-3t+3y 


2  Aj  -  f(xj)  +  £»(xj) 


-  d. 


whence 


0  i  9(t  -  5)A^  i 


Coeibinittg  this  with  (5.8).  we  get  that 


Q.BJ). 


Pfitmt  SX  3!!heere«  As  ateted  in  Section  4.  we  eowibine  e 

coapactneee  nxinwent  with  induction  to  proof  this  reeult. 

each  point  t  •  [a.b].  we  chooee  a  >  0  that  detenines  an  open 
intexTal  •  (t-b^.t+b^).  Since  (I^)  feme  an  open  cowering  of  the 
conpact  interral  [a.b].  there  existe  a  finite  enbeoTering  of  [a.b]. 
Horeorer.  for  h  -  aaxlh^)  aafficientlp  tBall,  each  interral 

one  of  the  interrale  of  the  enbeorering.  The  proof 
relies  hearily  190a  exploiting  the  local  properties  of  f  on  each  interral 
of  the  finite  enbeorering. 

In  choosing  b^,  we  consider  four  oases. 
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1.  If  f'(t)  0,  then  oKooae  6^  >  0  sneh  that 

0  <  f'(x)  <  3f*(y) 
for  all  X,  7  a  n  [a,b]. 

2.  If  f'(t)  *  0  bat  f"(t)  ^  0,  tbaa  obooae  6^  >  0  aoeb  that 

0  <  f*'(x)  <  1.5f»'(y) 
for  all  X,  7  a  I^  0 

3.  If  f'(t)  “  f**(t)  ■  0  bat  f^^Nt)  ^  0,  then  obooaa  >  0  aaeb  that 

0  <  f^^^x)  <  l.lf^^^7) 

for  all  X,  7  a  I^  ^  [t,b], 

(3i 

4.  If  f'(t)  “  f"(t)  •  f'  ^(t)  «•  0,  tbaa  obooae  B^  aoeb  that,  for  all 
7  a  [t,t^^)  n  [a,b],  aitber 

a.  f'(7)  ■  0  or 

b.  for  aoaa  eoaataata  aad  r» 

mj(r-t)*  i  r(7) 

•kara  ^  ^  ^  ^  S  ^ 

aad,  for  all  7  a  (t-B^,t]  H  U»b],  aitbar 

a.  f'(7)  -  0  or 

b.  for  aoaa  eoaataata  aad  a. 


■jCt-y)*  i  f’(7)  i«4(t-7)* 


where 

It  ie  poseible  to  ehoose  6^  to  sotiefj  Case*  1-3  beeenee  the  first  three 
derlTstlwes  of  f  are  eoatlnooiM.  If  f^^^(t)  ^  0,  tbes  Case  4  follow*  frow 
the  eoatinnity  of  f^^^.  Otherwise*  it  follows  directly  froa  Assaaptioa  2 
of  Theorea  4.1. 

To  prove  that  the  iadaotioa  hjrpothesis  (stated  ia  the  abbreviated 
proof  of  Theorea  4.1  ia  Seetioa  4)  reaaias  valid  whea  aa  approxiaate 
derivative  valae  d^  is  aodified  ia  Step  3  of  the  Ezteaded  Two-Sweep 
Algoritha,  we  eoasider  a  aoaber  of  eases  depeadiag  apoa  the  properties  of  f 
at  t*  where  ^  iaterval  aader  eoasideratioa.  We  prove 

the  last  ease  ia  the  iadaotioa  hypothesis  first. 

Case  1:  Asaaae  that  ^  f'(t)  4  0. 

Case  1.1:  Aasasie  that  u  u 

decreased  to  d^  oa  the  Forward  Sweep.  Heaee,  d^  2.  2.  Siaee 

A^_^  >  f*(y)  for  soae  y  s  it  follows  free:  the  choice  of  1^  that 

f'{Xj)  i  Therefore, 

0  i  d*  -  f*(Xj)  i  4^  -  f’(Xj)  i  oh*. 

Case  1.2:  Assoae  that  **4  4^  is  iacreas*4  to  4^  oa 

the  Porwar4  Sweep.  If  4^  i,  f'(x^),  thea 

0  i  f'(Xj)  -  4*  i  f'(Xj)  -  4j  i  oh*. 

(hi  the  other  haa4,  if  4^  2  f*(x^),  thea 
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0  i  dj  -  f'(x^)  i  d*. 

To  boBsd  d*,  aoto  tkat  i  liy  tho  eAoiee  of  I^.  Tkoxoforo. 

SAf  i  dj^^  i  SAj  ♦  ch^. 

2 

Tkii  iaoqulity  togotAor  vltA  tAe  oboorratloA  tbat  tbo  eorre  x  >  (y^S)  it 
eostftisod  lx  H  fox  3  i  y  i  d,  thamt  tAxt  d^^  i  (eh  )^»  tlxoo  •  Aj* 

Heoee*  for  h  txffieioxtly  mall* 

If(xj)  -  d*l  i  oh®. 

Cate  1.3:  Sixoe  **<*j_i)  i  ^A^-i  “*  f’(x^)  i  3Aj.  a  tiailax 
axgaaeBt  thovt  that 

If(Xj)  -  d*l  i  eh® 

after  the  Baekward  Sveep. 

Cate  2:  Attne  that  ^  f"(t)  d  0. 

(lx  thit  oate«  t  axat  be  oxe  of  the  exdpoixtt  of  the  ixterral  [a,b],  tlxee 
otherwlte  f  voxld  xot  be  aoxotoxe.)  At  lx  (^te  1.  the  choice  of  extxret 
that 

f'(Xj.j)  i  3Aj,j,  f»(Xj)  i  3Aj_j. 

f'(Xj)i3Aj,  r(Xj^j)iSA^. 

Therefore*  a  tiailar  argxaext  ahewt  that 

If'(Xj)  -  d*l  i  ch® 


•t  tbe  taxBinatiOB  of  Stop  3  ia  tkis  cost  ••  well 


Cose  3:  Assaae  tfcot  ^  ■  f'*(t)  •  0, 

#  0  sad  i  t. 

Cose  3.1:  Assaae  that  *  ^-1  ^  ^i-1  ^-1 

deereosed  to  d^.  Froa  tke  ehoioe  of  it  follows  tAot  f"(x)  <  0  for 
X  s  oad  X  <  t.  Therefore,  f*(x^)  Heaee,  os  ia  Cose  1.1, 

0  i  d^  -  f'{x^)  i  d^  -  f*(x^)  i  ok*. 

Cose  3.2:  Assaae  tkot  iaereosed  to  d^. 

Agoia,  siaee  f"(x)  <  0  for  x  s  oad  x  <  t,  it  follows  tkot 
f'(x^^^)  i,  A^.  Coaseqaeatlp,  tke  orgaaeat  ased  ia  Cose  1.2  skows  tkot,  for 
k  saffieieatly  saoll, 

lf»(Xj)  -  d*i  i  ok^ 

ia  tkis  ease  os  well. 

Cose  3.3:  Assaae  tkot  u  oad  d^  is  deereosed  to  d^. 

If  f»(Xj)  i  dj,  tkea 

0  i  d*  -  £'(Xj)  i  -  f»<x^)  i 

Tkerefore,  assaae  tkot  f'(x^)  >  d^,  oad  let  t  ”  S^xee 

A^  •  f  (x)  for  soae  x  s  (x£*X|4.2l  f*(x)  ^  l.lf’Cy)  for  oil  x,  y  s  X^, 

it  follows  froa  tke  Taylor  series  expoasioas  of  f*(x^)  oad  f*(x)  okoat  t 
tkot 

f*(Xj)/Aj  -  f**Nx)(Xj-t)*/f^*Ny)(r-t)*  i  1.1(t+1)*/t*. 
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CoM*4Matl7.  if  T  2.  vk«r«  t^Mi  f  (x^)  i  8A^  i  d^,  a 

eoatradietioa.  Tkmafora*  0  ^  Y  ^  l/(«r-l).  A  alaple  ealexlxtiea*  aiailar 
to  tAa  oaa  aaad  ia  tAa  proof  of  LoHaa  5.1,  aAova  that 

3Aj  i  d*  i  f(Xj)  i  SAj  +  dlf^^^l.h*. 

^aaeo 

0  i  f*(Xj)  “  d^  i 

la  aeeordaaee  vitk  tka  iadaetioa  kppotkaaia. 

Caaa  3.4:  Aaaaaa  tkat  a  **4  4^  ia  iaereaaed  to  d^. 

If  d*  i  tAaa 

0  i  f»(x^)  -  d*  i  f(x^)  -  d^  i  aaxie,  dlf^^^ljk®. 

Oa  tka  otkar  kaad,  if  d^  2.  fCa^),  tkaa  aa  argaaaat  aiailar  to  tka  oaa  aaad 
ia  Caaa  3.3  togatkar  witk  tka  iadaetioa  kypotkaaia  akowa  tkat 

Vi  i  *^i-i  * 

aw  9 

akara  e  ia  giaaa  ia  (4.2).  Tharafora,  aiaea  y  *  (z-3)  ia  ooataiaad  ia  K 
for  3  i  X  i  4  aad  ^4^^^,dj)  a  it  followa  tkat 

0  i  d*  -  f»<x^)  i  4j  i  (dlf^^^l.  ♦  c)V, 

akiak,  for  k  aaffioiaatly  aaall,  aatiafiaa  tka  iadaetioa  kypotkaaia. 

Caaa  4:  Aaaaaa  tkat  c  I^,  f'(t>  -  f»»(t)  •  0, 

f^^Nt)  d  0  aad  t  1  Aa  argaaaat  aiailar  to  tka  oaa  aaad  ia  Caaa  3 
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ahemt  that  tlie  iateetion  hjrpotltasi*  kolds  ia  this  ease  as  well. 

Case  5;  Asaaae  that  c  f'(t)  -  £"(t)  -  £**\t)  -  0 

aad  it,  la  this  ease,  either  f*(j)  ■  0  for  all  y  i  t  ia  or  f*(y) 
satisfies  the  boaad  ia  Coaditioa  4b  oa  If  f'(y)  ••  0,  thea  both 
aad  are  aero.  Heaee,  if  d^  ^  0,  thea  4^  is  set  to 

aero  oa  the  Forward  Sweep  of  the  Esteaded  Two-Sweep  Algoritha. 

Fartheraore,  siaee  d^  is  aot  modified  agaia,  d^  •  f'(x^)  »  0  at  the 
teraiaatioa  of  the  Step  3.  Therefore,  assaae  that  f'(y)  satisfies  the 
boaad  ia  Coaditioa  4b  oa  throaghoat  the  reaaiader  of  this  ease. 

Case  5.1:  Assaae  that  u  u  aad  d^  is 

decreased  to  d^.  Thea,  siaee  A^_j  «  f*(y)  for  some  y  s 
follows  froa  Coaditioa  4b  oa  that 

f*(Xj)/Aj^j  i  a^(t-Xj)*/aj(t-y)*  i  1,1. 

Therefore,  siaee  d^  2,  d^  2, 

©id*-  f'(Xj)  i  4j  -  f'(Xj)  i  eh^. 

Case  5.2:  Assaae  that  **4  ©£  f*  iaereased  to  d^.  Aa 

argaaeat  siailar  to  the  oae  abowe  shows  that  f'(x^^^)  i  I.IA^. 
Coaseqaeatly,  as  ia  Case  1.2, 

If'(Xj)  -  A*\  i  eh^ 

for  h  aaffieieatly  aaall. 


Case  5.3:  Assaae  that  •  &£  u  *>4  4^  is  deereased  to  d^ 


H*ao*.  if  f'(x^)  i  tkea 


0  i  dj  -  f'(Xj)  i  d^  “  f'(x^)  i  ck^. 

Ik*r*fore,  •••«■•  that  f  (z^)  >  d^*  lAeaea 

0  i  f(xj)  -  dj  i  £*(xj). 

To  boaad  £*(Xj),  lot  t  -  (t~x^^j)/hj.  Thoa,  aiaoe  -  £'(x)  for  aoae 
X  .  [x^.x^^jl, 

f'(x^)/Aj  i  a^(x^-t)*/aj(x-t)*  i  1.1(t+1)*/t* 

by  Coaditioa  4b.  Coaaoqaoatly,  if  y  i,  1/<m->1),  ahore  a*  ■  thoa 
f'(x^)  i,  3A^  i  d^,  a  ooatradiotioa.  Thozaforo.  0  £  y  ^  l/(«i-l).  Baaee,  if 
a  >  3,  thaa  f*(Xj)  «  o(h^)«  aad  the  iadaetioa  hypotheaia  holda  for  h 
aaffioieatly  aaall.  Oa  the  other  haad,  if  a  •  3,  thea  e^aadiag  f'(x^)  aa 
a  Taylor  aeriea  about  t  ahova  that 

f»(x^)  -  |f^^\x)(xj-t)^  i  ilf^^^lji*. 

aa  required. 

Caae  5.4:  Aaauaa  that  ^d^_j,dj^)  a  xxd  d^  ia  iaeraaaad  to  d^. 
Thea.  if  i*  i  £'<Xj), 

0  i  f'(Xj)  -  d*  i  f»(Xj)  -  d^  i  auxCe.  tlf^^^ljh*. 

to  the  other  head,  if  d*  i  £'(x^),  let  y  ■  (t-x^)/h^_j.  Thaa.  aa  arguaaat 
aijailar  to  the  oae  abore  together  with  the  iadaetioa  hypotheaia  ahowa  that 
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for  T  2.  vKore  «*  ■  Eomeo*  wo  ogoia  koro  tkot 

0  i  d*  -  i  d*  i  (2h^)*. 

Coaworoely,  if  0  1  y  1  l/(«i-l),  tkoa,  for  s  >  3» 

■  f'(x)  i  ■^(T+l)*k*  •  o(k^)* 

iddle,  for  t  *  3, 

•  f*(x)  -  |f^^Ny)(x-t)^  i  Slf^^^l^h*. 
la  oither  ease, 

0  i  d*  -  f»(x^)  a*  i  |a^_j  i  4lf^^^l,k^ 
for  k  oaffleioatly  wall. 

Caoo  (:  Aaaaao  tkat  c  f*<t)  •  f"(t)  -  f^*^t)  -  0 

aad  t  ^  ^  alallar  argaaoat  to  tko  oao  aaod  ia  Caoo  5  akowo  tkat  tko 

iadaetioa  kypotkooio  koldo  ia  tkio  oaao  aa  wall. 

Caoo  7:  Aaoaao  tkat  «  1^*  -  f**Nt)  -  0 

aaA  t  0  proof  of  tko  iadaetioa  kypotkooio  foliowx  oatily 

frOB  tko  obaorwatioa  tkat  f'(x^),  A^_,j  aad  A^  aro  oaek  boaadod  by 

lf<^>iy. 

Tkio  eoBplotoa  tko  proof  of  tko  tkird  oaao  of  tko  iadaetioa 
kypotkooio.  Wo  aow  ocaoider  tko  firat  two  oaaoo. 

Caoo  t:  Aooaao  tkat  ®  ^*<t)  ■  f"(t)  •  Q. 


I9\ 

V  ^(t)  ^  0  aad  t  t  First  aote  tlut*  for  li  szfficieatly  saall, 

Iz^_2**£]  ^  ^t*  **  slailsr  to  tks  oss  prosostod  is 

Csss  3*1  skows  that  1  ^£.2*  vkieb  it  follows  tliat*  if 

**i-2**i-l^  •  ^2  y  fii-2  U  ^-2* 

0  i  d*_2  “  i  ^i_l  “  ^*<*i-l^  i 

CoBsstasBtly,  d^^^  sstisfiss  (3.5)  at  tks  start  of  the  Forward  Swoop  for 

-i* 

Case  8.1.1:  Assoao  that  *  ^-1*  doeroaasd 

to  d^  only  if  d^^^  ^•*b  iaeroasod  to  d^_2  sad  oitbor 

1.  ^^i-2»^l-l^  boaadary  of  u  Ilj^_2  “  ^-2*  “ 

2.  boaadary  botwooa  *Bd  S^.^* 

la  tbo  first  easo,  d^^^  ^^£.2*  **  FB^^ioasly  aoatioaod. 

f*(*2_2)  i  Ai«2*  ^  ^i-*l*  ®^***^®B*»  ^  i<os«a  5.2, 

If(zj)  -  dji  iaazic,  jJlf^^^ljb*. 

Oa  tbo  otbor  bsad,  if  ®*  boaadary  botwooa  ^2.2  ^i-1* 

tboa  tbo  followiag  easo  applios  aftor  aotiag  that 
tbo  boaadary  of  tbsa  (^2.2*^i^  **** 

0 

Caso  8.1.2:  Assaao  tbst  *  ^1-1  ^  ^-2  ^i  ^®'*®**^ 

to  d^.  A  siaplo  ealealatioa  sbows  that  tbo  wortieal  distaaeo  froa 
(d^^^.d^)  to  tbo  boaadary  of  u  ^***  *8B*1  bo  2.75 

tiaos  tbo  aiaiaaa  distaaeo  froa  ^^2.2*^!^  ^>1* 

iaoqaalitios  (5.1),  (5.2)  aad  tbo  orror  boaads  oa  d^^^  *b^  ^2*  bt  follows 


-  4«  - 

tliat  tha  diataaea  froa  boaadary  of  la  laaa  tbaa  ox 

aqaal  to 

(2e  + 

Coaaaqaaatly* 

If(x^)  -  d*\  i  (6.5c  + 

Caaa  8.2:  Aaaaae  that  iaaraaaad  to  d^. 

If  d^  ±  f'(z^)«  thaa  we  agaia  hawe  that 

0  i  f*(Xj)  -  d*  i  f'(Xj)  -  dj  i  «.5a  +  |lf^^^ljh*. 

(hi  the  other  haad*  if  d^  },  f*(x^)*  thaa 

0  i  d^  -  f»(xj)  i  d^. 

Baeaaaa  t  a  (x^_,^,x^].  aa  azgaaeat  ciailar  to  the  oaa  praaaatad  la  Caaa  3.2 
ahowa  that 

*1+1  ^  ^'^*1+1^  +  eh^  i  3A^  ♦  6|f**^l,h®  +  eh* 
aad 

d*  i  (6lf<*^l.  +  a)V. 

whlah  ao^lataa  the  aaalyala  of  thla  oaaa. 

Caaa  S.S:  Aasaae  that  a  u  *>3  that  d^  la  doaraaaod 

to  d^.  Ihoxafoxa,  d^  2.  d^  1  SA^.  lowawax,  alaaa  t  a 


f'(x^)  i  A^.  Baaaa 
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0  i  -  f»(x^)  i  d^  -  f(*^)  i  («.5«  ♦  ||£<^^|Jk^ 

Case  8.4:  Assaaa  tkat  <dj_j,d^)  t  aad  tkat  d^  it  iaexaatad  to 
d^.  If  d*  i  f(x^).  tkaa 

0  i  f*(Xj)  -  4*  i  f(x^)  -  d^  i  «.5c  +  ||£^^>|Jk^ 

Ikaraforo.  aata.t  tkat  d^  I  f(x^).  la  additioa.  aott  tkat.  if 

«  Bi.j.  tkoB  «t  eoald  aot  kart  kad  (d^.^.d^)  t  4^,^  oa  tke 
Fonrard  Swaap.  Tkarafoxa.  tka  boaad 

•*1-1  "  i 

aataklitkad  at  tka  bagiaaiat  ot  Cata  8  atill  koldt.  MoraoTar.  tiaea  tka 

tlopa  of  tka  earra  tkat  foxat  tka  boBadarj  bataaaa  K  aad  £  it  latt  tkaa  or 
•«aal  to  oaa. 

Old*-  f»(x^)  i  (C  4  5J|f<^>|Jh^ 

Cata  9:  Attaae  tkat  c 

f  0  aad  t  t  Ix£»*£+i). 

Cata  9.1:  Attaaa  tkat  <dj_j.dj)  a  4^.^  U  u  ^ 
daoraatad  to  d*.  Ihoxafora.  d^  2  d*  i.  JA^.j.  Bowvwn,  f(x^)  i  A^.^  by 
tka  ekoiea  of  I^.  Haaoa. 

©id*-  f»(xj)  i  d^  -  f»(*j)  i  ok*. 

Cata  9.2:  Attaaa  tkat  <dj,d^^j)  a  4^  aad  d^  it  iaeraatad  to  dj.  If 
dj  i  f'(*j),  tkaa 

0  i  f'(*j)  -  d*  i  f'(Xj)  -  d^  i  ok*. 


Therefore,  aseeae  that  Note  that  the  isrerae  of  the  elope  of 

the  enrre  that  forma  the  hoaadarp  hetveea  H  amd  A  1*  !**•  thaa  or  equal  to 
’Oae.  Therefore,  as  ia  Casa  S.4, 

0  i  d*  -  f*(x^)  i  (c  ♦  jJlf^^^ljh*. 

Case  9.3:  Aasoae  that  u  E^.  Siaee  d^  i  3A^.  d^  eoald 

aot  haTe  beea  modified  ia  Caae  9.2.  laaee,  d^  mast  still  satisfy  (3.5). 
Coasider  the  folloviag  tvo  saheaaes. 

Case  9.3.1:  Assome  that  *  E^*  Note  that  d^  is  decreased 

to  d^  oaly  if  d^^^  was  iaereased  to  d^^^  aad  either 

^*i+l**i+2^  boBBdary  of  or 

2.  ^^£*^£42^  hoaadary  hetweea  aad  E^* 

la  the  first  ease,  d^^^  2.  ^£4X*  ^  additioa,  f*(x^^^)  ^  A^^^  by  the  choice 

of  I^.  Therefore,  f'(x^_^j)  i  d*^j  aad 

If(x^)  -  dji  ImaxCc, 

by  Lava  5.2.  Qa  the  other  haad,  if  is  oa  the  boaadary  betveea 

aad  E^*  thea  the  followiaq  ease  applies  after  aotiaq  that  (^£*^£42) 
closer  to  the  boaadary  of  thaa  (d£*4£4£)  *xs. 

Case  9.3.2:  Assame  that  ^^£»A£4£)  •  &£  oad  that  d^^  is  decreased  to 
d^.  As  ia  Case  8.1.2,  aote  that  the  horisoatal  distaaee  from  (^£*^£42)  ^o 
the  boaadary  of  H£  i*  !•••  thaa  or  eqaal  to  2.75  times  the  miaimw  distaaee 
from  (^£*^£42)  to  the  boaadary  of  E^.  Horeorer,  iaeqaalities  (5.1),  (5.2) 
aad  the  iadaetioa  hypothesis  oa  the  error  ia  d^^^^  imply  that  the  distaaee 
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froB  ^^£*^£4.2^  of  is  isss  than  or  o«ul  to 

(7.5c  + 

CoMoqoo&tly. 

If(*j)  -  dji  i  (fflf^^^l.  +  ^>k®* 

Cate  9.4:  Aacaac  tlut  (^£_£*^£)  *  ^>1  i^sraaaed  to  d^. 

If  d^  i,  f'(z^)*  tAaa  we  again  Aawe  that 

0  i  f*(x^)  -  d*  i  £*(Xj)  “  dj  i 

Qa  tAc  otAer  Aaad.  if  d^  1  f*(z^).  tAea 

0  i  d*  -  f»(xj)  id*  i  ♦  Bax{c.8)l£^^^|,})V, 

wAieli  followa  froa  aa  argaaeat  aiailar  to  the  oae  aaed  ia  Caae  3.4  after 
aotiag  tkat 

4£_j  i  3Aj_j  +  +  aax{c.«)lf^^^lj)h*. 


